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Both orbital and rotational dynamics employ the method of variation of parameters. We 
express, in a non-perturbed setting, the coordinates (Cartesian, in the orbital case, or Eule- 
rian in the rotation case) via the time and six adjustable constants called elements (orbital 
elements or rotational elements). If, under disturbance, we use this expression as ansatz 
and endow the "constants" with time dependence, then the perturbed velocity (Cartesian or 
angular) will consist of a partial derivative with respect to time and a so-called convective 
term, one that includes the time derivatives of the variable "constants." Out of sheer conve- 
nience, the so-called Lagrange constraint is often imposed. It nullifies the convective term 
and, thereby, guarantees that the functional dependence of the velocity upon the time and 
^ ■ "constants" stays, under perturbation, the same as it used to be in the undisturbed setting. 

The variable "constants" obeying this condition are called osculating elements. Otherwise, 
they are simply called orbital or rotational elements. 

When the dynamical equations, written in terms of the "constants," are demanded to 
be symplectic (and the "constants" make conjugated pairs Q, P), these "constants" are 
called Delaunay elements, in the orbital case, or Andoyer elements, in the rotational case. 
The Andoyer and Delaunay sets of elements share a feature not readily apparent: in certain 
cases, the standard equations render these elements non-osculating. 

In orbital mechanics, the elements, calculated via the standard planetary equations, come 
out non-osculating when perturbations depend on velocities. This complication often arises 
but seldom gets noticed. To keep elements osculating under such perturbations, extra terms 
must enter the equations, terms that will not be parts of the disturbing function (Efroimsky 
& Goldreich 2003, 2004). In the case of parametrisation through the Kepler elements, this will 
merely complicate the equations. In the case of Delaunay parametrisation, these extra terms 
will not only complicate the Delaunay equations, but will also destroy their canonicity. Under 
velocity-dependent disturbances, the osculation and canonicity conditions are incompatible. 

Similarly, in rotational dynamics, the Andoyer elements come out non-osculating when 
the perturbation depends upon the angular velocity of the top. Since a switch to a non- 
inertial frame is an angular-velocity-dependent perturbation, then amendment of the dy- 
namical equations by only adding extra terms to the Hamiltonian makes these equations 



1 



render non-osculating Andoyer elements. To make them osculating, extra terms must be 
added to the dynamical equations (and then these equations will no longer be symplectic). 

Calculations in terms of non-osculating variables are mathematically valid, but their 
physical interpretation is problematic. Non-osculating orbital elements parametrise instan- 
taneous conies not tangent to the orbit. Their inclination, the non-osculating i, may differ 
much from the physical inclination of the orbit, given by the osculating i. Similarly, in the 
case of rotation, non-osculating Andoyer variables do correctly describe a perturbed spin but 
lack simple physical meaning. The customary expressions for the spin-axis' orientation, in 
terms of the Andoyer elements, will no longer be valid, if the elements are non-osculating. 
These expressions, though, will stay valid for osculating elements, but then the (correct) 
dynamical equations for such elements will no longer be canonical - circumstance ignored in 
the Kinoshita-Souchay (KS) theory which tacitly employs non-osculating variables. While 
the loss of osculation will not influence the predictions for the figure axis of the planet, it 
considerably effects the predictions for the orientation of the instantaneous axis of rotation. 

1 Kepler and Euler 

In orbital dynamics, a Keplerian conic, emerging as an undisturbed two-body orbit, is regarded as 
a sort of "elementary motion," so that all the other available motions are conveniently considered 
as distortions of such conies, distortions implemented through endowing the orbital constants 
Cj with their own time dependence. Points of the orbit can be contributed by the "elementary 
curves" either in a non-osculating fashion, as in Fig. 1, or in the osculating way, as in Fig. 2. 

The disturbances, causing the evolution of the motion from one instantaneous conic to another, 
are the primary's oblateness, the gravitational pull of other bodies, the atmospheric and radiation- 
caused drag, and the non-inertiality of the reference system. 

Similarly, in rotational dynamics, a complex spin can be presented as a sequence of configu- 
rations borrowed from a family of some elementary rotations. The easiest possibility here will be 
to employ in this role the Eulerian cones, i.e., the loci of the rotational axis, corresponding to 
non-perturbed spin states. These are the simple motions exhibited by an undeformable free top 
with no torques acting thereupon. 1 Then, to implement a perturbed motion, we shall have to go 
from one Eulerian cone to another, just as in Fig. 1 and 2 we go from one Keplerian ellipse to 
another. Hence, similar to those pictures, a smooth "walk" over the instantaneous Eulerian cones 
may be osculating or non-osculating. 

The physical torques, the actual triaxiality of the top, and the non-inertial nature of the refer- 
ence frame will then be regarded as perturbations causing the "walk." The latter two perturbations 
depend not only upon the rotator's orientation but also upon its angular velocity. 

2 Delaunay and Andoyer 

In orbital dynamics, we can express the Lagrangian of the reduced two-body problem via the 
spherical coordinates qj = {r , ip , 9} , then calculate their conjugated momenta pj and the 
Hamiltonian H(q, p) , and then carry out the Hamilton- Jacobi procedure (Plummer 1918), to 
arrive to the Delaunay variables 

{Qi, Q2, Q 3 ; Pi, P 2, P 3 } = {L,G,H;l,g,h} = 

1 Here one opportunity will be to employ in the role of "elementary" motions the non-circular Eulerian cones 
described by the actual triaxial top, when this top is unforced. Another opportunity will be to use, as "elementary" 
motions, the circular Eulerian cones described by a dynamically symmetrical top (and to treat its actual triaxiality 
as another perturbation). The main result of our paper will be invariant under this choice. 
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(1) 

{ y/JIa , J fia(l — e 2 ) , y/ia (1 — e 2 ) cos « ; — M Q , — , — } , 

where /i denotes the reduced mass. 

Similarly, in rotational dynamics one can define a spin state of a top by means of the three 
Euler angles qj = ip , 6 , <p and their canonical momenta pj , and then perform a canonical trans- 
formation to the Andoyer elements L,G,H,l,g,h.A minor technicality is that, historically, 
these variables were introduced by Andoyer (1923) in a manner slightly different from the set of 
canonical constants: while, for a free rotator, the three Andoyer variables G , H , h are constants, 
the other three, L , I , g do evolve in time (for the Andoyer Hamiltonian of a free top is not zero, 




Fig.l. The perturbed trajectory is a set of points belonging to a sequence of 
confocal instantaneous ellipses. The ellipses are not supposed to be tangent, nor 
even coplanar to the orbit at the intersection point. As a result, the physical velocity 
r (tangent to the trajectory at each of its points) differs from the Keplerian velocity 
g (tangent to the ellipse). To parametrise the depicted sequence of non-osculating 
ellipses, and to single it out of all the other such sequences, it is suitable to employ 
the difference between r and g , expressed as a function of time and six (non- 
osculating) orbital elements: 

${t, d , C 6 ) = f(t, Ci , C fl ) - g(t, Ci , C 6 ) . 
Since 

r ~ dt + 4-i dt C J - g + 4-, Qt c > > 

J=l J=l 

then the difference $ is simply the convective term (dr/dCj) Cj which 
emerges whenever the instantaneous ellipses are being gradually altered by the per- 
turbation (and the orbital elements become time-dependent). In the literature, 
<&(i , C\, . . . , Cq) is called the gauge function or gauge velocity or, simply, gauge. 




Fig.2. The perturbed trajectory is represented through a sequence of confocal 
instantaneous ellipses which are tangent to the trajectory at the intersection points, 
i.e., are osculating. Now, the physical velocity r (which is tangent to the trajectory) 
will coincide with the Keplerian velocity g (which is tangent to the ellipse), so that 
their difference <&(t C\ , . . . , Cq) vanishes everywhere: 

6 BC 

Ci C 6 ) = f(t, d , C 6 ) - g(t, d, C 6 ) = "£-Qf C i = 

This equality, called Lagrange constraint or Lagrange gauge, is the necessary and 
sufficient condition of osculation. 
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but a function of I , L and G). This way, to make our analogy complete, we may carry out 
one more canonical transformation, from the Andoyer variables {L , G , H , / , g , h} to "almost 
Andoyer" variables {L Q , G , H , l Q , g a , h} , where L a , l and g D are the initial values of L , I 
and g . The latter set consists only of the constants of integration; the corresponding Hamiltonian 
becomes nil. Therefore, these constants are the true analogues of the Delaunay variables (while 
the conventional Andoyer set is analogous to the Delaunay set with M used instead of M Q .). The 
main result obtained below for the modified Andoyer set {L , G , H , l a , g Q , h} can be easily 
modified for the regular Andoyer set of variables {L , G , H , I , g , h} . (Efroimsky 2005b) 
To summarise this section, in both cases we start out with 

dH {o) . _ dH® 

^ dp ' dq 

q and p being the coordinates and their conjugated momenta, in the orbital case, or the Euler 
angles and their momenta, in the rotation case. Then we switch, via a canonical transformation 



q = f(Q, P, t) 

p = X (Q, P, t) 



(3) 



to 



ffu* FFH* 

where Q and P denote the set of Delaunay elements, in the orbital case, or the (modified, as 
explained above) Andoyer set { L Q , G , H , l Q , g Q , h } , in the case of rigid-body rotation. 

This scheme relies on the fact that, for an unperturbed Keplerian orbit (and, similarly, for an 
undisturbed Eulerian cone) its six-constant parametrisation may be chosen so that: 

1. the parameters are constants and, at the same time, are canonical variables { Q , P } with a 
zero Hamiltonian: 7i*(Q, P) = ; 

2. for constant Q and P , the transformation equations are mathematically equivalent to 
the dynamical equations (j2J). 

3 When do the elements come out non-osculating? 
3.1 General-type motion 

Under perturbation, the "constants" Q , P begin to evolve so that, after their substitution into 

q = f{Q{t), P(t), t) 

(5) 

p = x (Q(t), Pit), t) 

(/ and x being the same functions as in © ), the resulting motion obeys the disturbed equations 



d(n {o) + ah) d(n {o) + ah) 

q = -± L , V = ~ ~ L ■ (6) 

op aq 

We also want our "constants" Q and P to remain canonical and to obey 

■ d{H* + AH*) ■ d(H* + AH*) 

Q ~ OP ' dQ { ] 
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where 



H* = and AH* (Q , P t) = AH{q{Q,P,t), p(Q,P,t), t) . (8) 

Above all, an optimist will expect that the perturbed "constants" Cj = Q\ , Q2 , Q3 , Pi , P% , P3 
(the Delaunay elements, in the orbital case, or the modified Andoyer elements, in the rotation 
case) will remain osculating. This means that the perturbed velocity will be expressed by the 
same function of Cj(t) and t as the unperturbed one used to. Let us check to what extent this 
optimism is justified. The perturbed velocity reads 

q = 9 + $ (9) 

where 

9(c(t) , t) s *2gM> (10) 

is the functional expression for the unperturbed velocity; and 

3=1 OU 3 

is the convective term. Since we chose the "constants" Cj to make canonical pairs (Q, P) obeying 
(0-IB1), with vanishing 7i* , then insertion of (J7J) into will result in 

•■£4^ + s^-^- (12) 

So the canonicity demand is incompatible with osculation. In other words, whenever a momentum- 
dependent perturbation is present, we still can use the ansatz (J3J) for calculation of the coordinates 
and momenta, but can no longer use (|T3jl for calculating the velocities. Instead, we must use (f*T3|) . 
Application of this machinery to the case of orbital motion is depicted on Fig.l. Here the constants 
Cj = (Q n , Pn) parametrise instantaneous ellipses which, for nonzero $ , are not tangent to the 
trajectory. (For more details see Efroimsky & Goldreich (2003).) In the case of orbital motion, 
the situation will be similar, except that, instead of the instantaneous Keplerian conies, one will 
deal with instantaneous Eulerian cones (i.e., with the loci of the rotational axis, corresponding to 
non-perturbed spin states). 

3.2 Orbital motion 

In the orbital-motion case, osculation means the following. Let the unperturbed position be given, 
in some fixed Cartesian frame, by vector function / : 

r = f(C u C 6 , t) , r = {x, y, z} . (13) 

Employing this functional ansatz also under disturbance, we get the perturbed velocity as 

F = g(C x (t) , . . . , C 6 (t) ,t) + $ (C x (t) , . . . , C 6 (t) , t) (14) 

where 

g-f and * s gg<7,. (15) 
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The osculation condition is a convenient (but totally arbitrary!) demand that the perturbed 
velocity r has the same functional dependence upon t and Cj as the unperturbed velocity g : 

f(d(t) , . . • , C e (t) , t) = /(Ci(f) , . . . , C 6 (t) , t) , 

(16) 

?(C x (t) , . . . , C 6 (t) , t) = g (d(t) , . . . , C 6 (t) , t) . 
or, equivalently, that the so-called Lagrange constraint is satisfied: 

oo o 

£ T^rQ = CeC*) , «) where * = . (17) 

Fulfilment of these expectations, however, should in no way be taken for granted, because the 
Lagrange constraint (|T7jl and the canonicity demand ([3- EJ) are now two independent conditions 
whose compatibility is not guaranteed. As shown in Efroimsky (2002a,b), this problem has gauge 

— * 

freedom, which means that any arbitrary choice of the gauge function <fr {C\{t) , . . . , C§(t) , t) 
will render, after substitution into 1)131- I14jl . the same values for r and r as were rendered by 
Lagrange's choice ()17j) . 2 As can be seen from (II 2jl . the assumption, that the "constants" Q 
and P are canonical, fixes the non-Lagrange gauge 

£ j^Cj = ^(Ci(t), Ca(f) , t) where $ = — £ . (18) 
j=i OL, j °P 

It is easy to show (Efroimsky & Goldreich 2003; Efroimsky 2005b) that this same non-Lagrange 
gauge simultaneously guarantees fulfilment of the momentum-osculation condition: 

r , . . . , C 6 (t) , t) = /(Ci(t) , . . . , C 6 (f) , t) , 

(19) 

p(d(t) , . . . , C 6 (t) , t) = g(d(t) , . . . , C 6 (t) , t) . 

Any gauge different from (|18jl. will prohibit the canonicity of the elements. In particular, for 
momentum-dependent A7i , the choice of osculation condition $ = would violate canonicity. 

For example, an attempt of a Hamiltonian description of orbits about a precessing oblate 
primary will bring up the following predicament. On the one hand, it is most natural and conve- 
nient to define the Delaunay elements in a co-precessing (equatorial) coordinate system. On the 
other hand, these elements will not be osculating in the frame wherein they were introduced, and 
therefore their physical interpretation will be difficult, if at all possible. Indeed, instantaneous 
ellipses on Fig.l may cross the trajectory at whatever angles (and may be even perpendicular 
thereto). Thence, their orbital elements will not describe the real orientation or shape of the 
physical trajectory (Efroimsky & Goldreich 2004; Efroimsky 2005a). 

For the first time, non-osculating elements obeying (f2*2"j) implicitly emerged in (Goldreich 1965) 
and then in Brumberg et al (1971), though their exact definition in terms of gauge freedom was 
not yet known at that time. Both authors noticed that these elements were not osculating. 
Brumberg (1992) called them "contact elements." The osculating and contact variables coincide 
when the disturbance is velocity- independent. Otherwise, they differ already in the first order of 
the velocity-dependent perturbation. Luckily, in some situations their secular parts differ only in 
the second order (Efroimsky 2005a), a fortunate circumstance anticipated yet by Goldreich (1965). 

2 Physically, this simply means, r on Fig.l can be decomposed into g and $ in a continuous variety of ways. 
Mathematically, this freedom reflects a more general construction that emerges in the ODE theory. (Newman & 
Efroimsky 2003) 
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3.3 Rotational motion 



where 



In rotational dynamics, the situation of an axially symmetric unsupported top at each instant of 
time is fully defined by the three Euler angles q n = 6 , , if) and their derivatives q n = 6 , <f> , if>. 
The time dependence of these six quantities can be calculated from three dynamical equations of 
the second order and will, therefore, depend upon the time and six integration constants: 

<?n = fn (S\ , . . . , Sq , t) , 

(20) 

Qn = g n (Si , . . . , S 6 , t) , 

the functions g n and /„ being interconnected via g„ = df n /dt , for n = 1,2,3 = if) , 9 , <fi. 
Under disturbance, the motion will be altered: 

q n = f n (Si(t), S 6 (t) , t) , 

(21) 

q n = g n (Si(t) , . . . , S 6 (t) ,*)+$„ (Si(t) , . . . , S 6 (t) , t) , 

6 df ■ 

^ n (S 1 (t),...,S e (t),t) = E^f^ • ( 22 ) 

Now choose the "constants" Sj to make canonical pairs (Q, P) obeying (JZI-EJ), with H* being 
zero for (Q, P) = (L , G , H , l , g , h) . Then insertion of (JJJ) into (j22J) will result in 

,. m) ,..., m , t)s ^ + E ^p = ?^lL , (23) 

so that the canonicity demand (j3-|HJ) violates the gauge freedom in a non-Lagrange fashion. This 
is merely a particular case of ()12j) . 

This yields two consequences. One is that, in the canonical formalism, calculation of the 
angular velocities via the elements must be performed not through the second equation of (|2Ujl 
but through the second equation of (}2T|) . with (J23|) substituted therein. This means, for example, 
that in Kinoshita (1977) expressions (2.6) and (6.24 - 6.27) render not the angular velocity of 
the Earth relative to the precessing frame (wherein the Andoyer variables were defined) but 
the angular velocity relative to an inertial frame. (For an extended explanation of this fact see 
Efroimsky 2005b.) This, however, is not a drawback of the Kinoshita-Souchay theory but rather 
its advantage, because it is the angular velocity relative to an inertial frame that is directly 
measurable. (Schreiber et al 2004) 

The second consequence is that, if we wish to make our Andoyer variables osculating (so that 
the second equation of (|2*U|) could be used), the price to be payed for this repair will be the loss of 
canonicity. (Angular-velocity-dependent disturbances cannot be accounted for by merely amend- 
ing the Hamiltonian!) The osculating elements will obey non-canonical dynamical equations. 

To draw to a close, we would add that, under some special circumstances, the secular parts 
of contact elements may coincide in the first order with those of their osculating counterparts. 3 
Whether this will be the case for the Earth or Mars remains to be investigated. This matter 
will be crucial for examining the validity of the presently available computations of the history of 
Mars' obliquity. 4 

3 In regard to orbital motions, this possibility was anticipated yet in 1965 by Peter Goldreich. As demonstrated 
by Efroimsky (2005a), this is true for constant rate of frame precession (but not for variable precession). 

4 The pioneer study on this topic was conducted by Ward (1973) in a direct manner and was, therefore, exempt 
from the problems associated with the loss of osculation. However, some of his successors chose the canonical 
formalism and exploited the Hamiltonian borrowed from the Kinoshita theory. 
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4 Conclusions 



In this talk we have explained why the Hamiltonian theory of Earth rotation renders non-osculating 
Andoyer elements. We have also explained how this defect of the theory should be mended. 

In attitude mechanics, osculation loss has the same consequences as in the theory of orbits: 
while this defect of the theory has no influence upon the theory's predictions for the figure axis 
of the planet, it considerably effects the predictions for the orientation of the instantaneous axis 
of rotation. 

In our paper Efroimsky (2005b) we shall demonstrate that, even though the Andoyer variables 
in the Kinoshita-Souchay theory are introduced in the precessing frame of the Earth orbit, they 
return the angular velocity not relative to that frame, but relative to an inertial one. This is not 
a drawback of this theory but rather its advantage, because it is the angular velocity relative to 
an inertial frame that is directly measurable at present. (Schreiber et al 2004) 
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